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We develop a technique for generating multi-photon nonclassical states via interference between 
coherent and Fock states using quantum catalysis. By modulating the coherent field strength, 
the number of catalyst photons and the ratio of the beam splitter upon which they interfere, a 
wide range of nonclassical phenomena can be created, including squeezing of up to 1.25 dB, anti- 
and super-bunched photon statistics and states exhibiting over 90% fidelity to displaced coherent 
superposition states. We perform quantum catalysis experimentally, showing tunability into the 
nonclassical regime. Our protocol is not limited by weak nonlinearities that underlie most known 
strategies of preparing multi-photon nonclassical states. Successive iterations of this protocol can 
lead to direct control over the weights of higher-order terms in the Fock basis, paving the way towards 
conditional preparation of "designer" multi-photon states for applications in quantum computation, 
communication and metrology. 

PACS numbers: 42.50 Ar, 42.50 Dv 



INTRODUCTION 

Quantum technologies promise enhanced performance 
of computation, communication, sensing, and simulation 
protocols. Photonic quantum information processing has 
been immensely successful in harnessing these advantages 
in several settings [ - ]. This is because the quantum 
states of photons can be prepared, processed and mea- 
sured with ease and precision whilst simultaneously pre- 
serving fragile quantum phenomena, even in a hostile en- 
vironment as they do not suffer from any significant cou- 
pling with the external environment. Depending on the 
optical degrees of freedom in which quantum informa- 
tion is stored, manipulated and detected, photonic quan- 
tum information processing typically operates in two dif- 
ferent regimes: discrete variable (DV) and continuous 
variable [4-6] (CV). In contrast to DV states residing in 
finite-dimensional Hilbert spaces, the CV quadratures of 
an optical field allow the encoding of information over 
infinite dimensions. This approach has already demon- 
strated advantages in the implementation of secure key 
distribution [ ], and there are ongoing efforts to exploit 
its potential in computation and communication [8-10]. 
For a protocol to provide quantum advantages not at- 
tainable classically, it is crucial that the quantum states 
involved exhibit non- Gaussian features in their Wigner 
distributions in phase space [ 1—13] . 

Generating non-Gaussian states normally requires 
nonlinearities of the third order in the field operators [4] . 
However, such effects are negligible at low photon flux. In 
this regime — the most relevant scenario in current exper- 
iments at the quantum level — a more feasible strategy in- 
volves the use of conditional probabilistic operations [14]. 
These comprise the manipulation of the quantum state 



by means of linear optical elements [ 5] , which gener- 
ates an effective nonlinearity by accepting only particu- 
lar outcomes of measurements on ancillary modes. Such 
schemes have allowed the production of superpositions of 
coherent states [ , ], and manipulation at the single- 
photon level [18-22]. 

In this Article, we investigate the possibilities offered 
to multi-photon quantum state engineering by a particu- 
lar interferometric scheme called quantum catalysis [14], 
which uses a photonic Fock state to modulate the prob- 
ability amplitudes of a coherent state. The analogy to 
catalysis is motivated by conditioning on the same num- 
ber of photons, say /c, as is initially interfered with the 
coherent state, as shown in Fig. 1. We generalise the 
result in Ref. [ ] (in which k = 1, a <C 1 and beam 
splitter transmissivity t <C 1) across a wide range of 
parameters and to higher-order in the photon number 
distribution, introducing a broad family of multi-photon 
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FIG. 1: A Photon Catalysed Optical Coherent state |PCOC) 
is produced following the interference between a coherent 
state \a) and a /c-photon Fock state \k) (at a beam split- 
ter of reflectivity |r| 2 = 1 — |t| 2 ), conditional on measuring k 
photons at one output port. 
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FIG. 2: Squeezing behaviour as a function of interaction parameters |r| 2 and a. The range of parameter space that generates 
squeezing is shown in (a), with the probability of successful catalysis shown in (b). The white lines show the parameters which 
generate the maximum squeezing of 1.25 dB below the shot-noise limit and the blue lines show the maximum non-Gaussianity, 
corresponding to a deviation from a minumum uncertaity state of 4.77 dB. (c) shows the quadrature variances, relative to the 
vacuum (unsqueezed) state in units of dB, at a coherent state amplitude of a = 1. Red dotted line, blue dashed line: X and 
Y quadrature variance, respectively. Green solid line: joint variance. 



quantum states. We experimentally measure the second- 
order autocorrelation function (0) for states catal- 
ysed by one photon (k = 1), as well as their photon- 
number distributions. The states we produce are deep in 
the quantum regime, and our strategy allows for the pro- 
duction of other nonclassical states such as squeezed and 
Schrodinger kitten states. In spite of injecting and ex- 
tracting the same number of photons, the photon- number 
distribution of the final quantum state is radically mod- 
ified due to quantum interference. This operates at two 
levels: firstly, each Fock state contributing to the coher- 
ent input state interferes with the catalysing k photons, 
undergoing a generalised Hong-Ou-Mandel type interfer- 
ence [ ]. Secondly, the outcomes of these individual in- 
terferences then add coherently due to the inherent phase 
reference in a coherent state. This combination results in 
a completely novel class of multi-photon quantum states 
with a whole gamut of uniquely nonclassical properties. 

INTRODUCING PCOC STATES 

The schematic for the generation of a photon-catalysed 
optical coherent (PCOC) state is shown in Fig. (1). The 
conditioned state at the output is given by 

co n 

|PCOC) = M a ,r Y\ ^r=C n (r,t,k) \n) , (1) 
c — l \/n\ 



where 
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r, t are the reflectivity and transmissivity of the beam 
splitter, respectively (satisfying |r| 2 + \t\ 2 = 1) and J\f a , r 
is the normalisation constant. The states cover a range 
going from a coherent state (r = 0) to a /c-photon Fock 
state (r = 1). By tuning the parameters of the inter- 
action, namely the beam splitter reflectivity r, coherent 
state amplitude a and the number of catalyst photons 
the coefficients of the photon-number terms (Fock layers) 
may be modulated, generating a wide range of nonclas- 
sical phenomena, as we show next. 

Squeezed states form an essential part of the quantum 
toolbox, and PCOC states can, for a range of parameters, 
yield an X quadrature variance A 2 X below the standard 
limit. This is shown in the dark regions of Fig. 2(a). 
Analytical expressions are obtained for the variance of 
the quadratures: 
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FIG. 3: Wigner functions of the photon catalysed state (a) and displaced coherent state superposition (CSS) (b), of amplitude 
a — 0.9, displaced by /3 = 0.8 from Eq. (8). Their Wigner functions are similar, indeed the photon-catalysed state exhibits a 
fidelity of 0.9 with respect to the displaced coherent state superposition. This result shows that CSS states can be generated 
without squeezing using the catalysis scheme. 




(a) k = 2, |r| 2 = 1/2, a = 2 (b) k = 3, |r| 2 = 1/4, a = 2 (c) k = 6, |r| 2 = 4/9, a = 2.7 



FIG. 4: Wigner functions of higher-order catalysed states showing the rich structure that can be achieved. 
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These quadrature variances, relative to their vacuum val- 
ues of 1/4, are plotted on a logarithmic scale (units of dB) 
in Figs. 2(a) and 2(c) (for a = 1) showing the range of 
squeezing for a and |r| 2 . 

By minimizing the expressions above for A 2 X, the 
largest squeezing is attainable is 3/16, below the vac- 
uum noise level of 1/4 by 1.25 dB. The maximum vari- 
ance of 3/4 corresponds to 4.77 dB antisqueezing. The 
parameters that yield these values are found by solving 
A 2 X=3/16 under the constraints < r < 1, \a\ <0 
(using Mathematica) . The solutions to this equation are 
not uniquely valued: they form two continuous lines in 



(|r| 2 ,a) parameter space, as shown in Fig. 2(a). The 
equations corresponding to the loci of minima are 
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Similarly for the locus of maximum non-Gaussianity is 
described by 



(6) 



Reaching the maximum squeezing and non-Gaussianity 
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by this scheme is therefore independent of the initial co- 
herent state amplitude a. 

However, the amount of squeezing obtainable, partic- 
ularly at high a, is offset by the probability of successful 
catalysis: if the coherent state is large, the probability 
of detecting one photon is small. The probability of suc- 
cessful catalysis is given by 



ftcoc O, r) = ^2 \ c "> ( a > r ) I 



n=0 

=e -r 2 \*\ 2 (! _ r 2 (! _ | a |2 ( r 2 (3 + | a |2J _ 2 _ r 4| Q |2))) 

(7) 

and it is shown in Fig. 2(b). For low a, catalysis is fairly 
likely. Indeed, at a = 1 and |r| 2 = 0.332, which gen- 
erates maximum squeezing, the probability of successful 
catalysis is about 47%. 

The behavior of both quadrature variances at a = 1 is 
shown in Fig. 2(c). The P quadrature is never squeezed, 
and the combined variances show that the state clearly 
satisfies the minimum variance condition. For |r| 2 < 0.2, 
the state asymptotically approaches a minimum uncer- 
tainty state, whereas this is lost for |r| 2 > 0.2. Since 
the state is pure, this means states in this region are 
non-Gaussian, implying negativity in their Wigner dis- 
tributions [23]. 

Non-Gaussian multi-photon states have been shown 
to be useful for computation [13, 24] and communica- 
tion [ ]. One class of non-Gaussian states of particular 
interest are coherent state superposition (CSS) states. 
As a superposition of coherent state amplitudes sepa- 
rated in phase space, they are the optical analog of the 
Schrodinger cat "dead" and "alive" state [ ] . When the 
separation between the amplitudes is small, the states 
are termed Schrodinger kitten states [ ] , and it has been 
shown that combining small CSS states can yield much 
larger CSS states necessary for universal quantum com- 
putation [28]. 

A CSS state of amplitude a, displaced in phase space 
by D (/}), is given by 



|CSS)=A/-£ 
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n=0 
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An example of a weak CSS state (a = 0.9), displaced in 
phase space (/? = 0.8), is shown in 3(b). With a partic- 
ular choice of interaction parameters (|r| 2 = 0.77, a = 
1.35), the PCOC state approximates a displaced weak 
coherent superposition state, with > 90% fidelity, with a 
Wigner function as shown in Fig. 3(a), which is clearly 
non-Gaussian. Previous schemes have typically required 
some initial squeezing [27, 29-32]; our results indicate a 
route to these states without such requirements. 

Simple extensions of the catalysis scheme allow for the 
preparation of more sophisticated quantum states. A cat- 
alyst of k > 2 photons in Eq. (1) leads to more elaborate 



Wigner functions as shown in Fig. (4). A further exten- 
sion is to concatenate iterations of catalysis. Each iter- 
ation i introduces an additional reflectivity n, such that 
the total reflectivity r = (7*1, . . . , ri). The state following 
/ iterations may be written 
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Such an iterative scheme forms the basis of a num- 
ber of conditional approaches towards quantum com- 
putation [ ]. The coefficients of the Fock layers can 
be directly controlled by the reflectivities at each itera- 
tion [ ]. As such, arbitrary quantum states can be gen- 
erated by heralding the successful catalysis interaction 
on a single mode, without invoking entanglement [34]. 
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FIG. 5: Experimental setup for probing single-photon catal- 
ysis. A Ti: Sapphire laser is doubled by a second harmonic 
generation crystal (SHG). The up-converted beam undergoes 
type-II co- linear parametric down- conversion (PDC), produc- 
ing pairs of orthogonally polarized photons in pure spectral- 
temporal states [35]. HWP1 ensures that the i7-polarized 
photon is detected by an avalanche photodiode (APD), which 
heralds the presence of a photon in the V polarization mode. 
The residual fundamental, which is used as the coherent state, 
is reflected by the dichroic mirror DM. This beam is then at- 
tenuated by a half- wave plate HWP2 and polarizing beam 
splitter PBS 2 to single-photon level intensity. The two beams 
are spatially recombined on PBS3, and synchronized by ad- 
justing the position of a delay stage. This operation is fol- 
lowed by interference on the variable beam splitter constituted 
by HWP3 and PBS4. The two polarization modes are then fil- 
tered by 3nm interference filters (IF), coupled to single-mode 
fibres and sent to photon-number-resolving time-multiplexed 
detectors (TMDs) for photon counting [36]. 



EXPERIMENTAL TUNING OF THE 
SECOND-ORDER COHERENCE 

To demonstrate the practical efficacy of our proposal, 
we experimentally generate PCOC states catalysed by 
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one photon. Using the setup shown in Fig. (5), we mea- 
sure the photon number distribution of the PCOC state 
at a range of beam splitter reflectivities. The variable 
beam splitter comprises a PBS (PBS3 - to spatially over- 
lap the coherent state and Fock states), a half- wave plate 
(HWP3) and interference beam splitter (PBS4). The 
beam splitter reflectivity and half- wave plate angle in the 
variable beam splitter are related by r = 1 + cos(20). 

The resulting modes after interference are sent to the 
detector array. By multiplexing APDs spatially and tem- 
porally, we are able to achieve photon number resolution 
up to eight photons plus vacuum [36]. Typical results 
from number-resolved counting are shown in Fig. (6), 
taking into account the binning of our detector. The 
action of catalysis as a Fock state filter appears in the 
modulation of the higher-order photon terms as the re- 
flectivity is tuned [21]. The single photon component 
is suppressed with respect to what one would expect in 
a coherent state. A similar effect is also present in the 
two- and three-photon components, albeit with a lower 
contrast. 

We also measure the second order autocorrelation func- 
tion (0) and show its variation with beam splitter re- 
flectivity |r| 2 . Fig. (7) shows that states lying below the 
classical bound (g^ 2 \0) = 1) are observed for high reflec- 
tivities. In this region, our experiment realizes a displace- 
ment of a single photon, as first observed in [ ]. This 
signature of nonclassicality disappears rather quickly as 
|r| 2 decreases below 0.9. At the other end where |r| 2 ~ 0, 
we clearly find coherence properties close to the ones of 
the input state \a). A significant feature in the curve is 
the presence of a peak centred around |r| 2 = 0.6, which 
attains a measured value of g^ 2 \0) = 2.5±0.23, above the 
limit of thermal states, which has been termed "super- 
bunching" [ ] . While this is not a signature of nonclassi- 
cality per se (in principle any value can be achieved by a 
mixture of coherent states), we can still make inferences 
of nonclassical features by observing the corresponding 
statistics measured in Fig. 6(a). In the ideal case, the 
Fock state filter should significantly suppress the single 
photon component for a reflectivity |r| 2 = \t\ 2 = 1/2, by 
Hong-Ou-Mandel interference [ ]. Its action on a coher- 
ent state should then bring it closer to a weakly squeezed 
state - which only contains contributions from even pho- 
ton numbers, and can reach #( 2 )(0)>3. While our de- 
tection scheme is not able to show squeezing explicitly, 
these photon counting measurements shows evidence of 
nonclassical statistics. 

Conclusions - We have introduced a general scheme 
for producing a broad class of useful multi-photon states 
without strong optical nonlinarities, whose generation is 
feasible with current technology. We have exhibited how 
a broad range of sophisticated quantum states can be de- 
signed starting from just single photons, coherent states, 
beam splitters and photon counters. In conjunction with 
quantum memories that can drastically increase multi- 




(a) Experimental data 




(b) Theoretical reconstruction 



FIG. 6: Measured (a) and theoretical (b) photon number dis- 
tributions as a function of beam splitter reflectivity |r| 2 for 
a coherent state mean photon number of \a\ 2 = 1.11. White 
shows the vacuum component, while red, blue and green rep- 
resent the one-, two- and three-photon terms, respectively. 
Note that the higher-order terms have been scaled for clar- 
ity; the captions X x Pij corresponds to a scaling factor X 
multiplied by the probability P of a click corresponding to i 
photons detected in TMD1 and j photons detected in TMD2. 



photon rates, and also act as tunable beam splitters [38], 
photon catalysis can be used as a building block in the 
construction of extended quantum networks [10]. 
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FIG. 7: Second order autocorrelation function g^ 2 \0) as a 
function of beam splitter reflectivity, for a coherent state mean 
photon number of \a\ 2 — 1.11. The points are measured val- 
ues based on the photon-statistics in Fig. 6(a), while the line 
is the expected value, incorporating inefficiency in our de- 
tection system. For \r\ 2 = 0,1, g^ 2 \0) approaches the val- 
ues for a coherent state and single photon, respectively. The 
regions showing anti-bunching g^ 2 \0) < 1 (strictly nonclassi- 
cal), bunching 1 < (0) < 2 and super-bunching g^ (0) > 2 
are clear. 
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